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Abstract 

We study the Tevatron search potential for minimal supergravity (mSUGRA) and find 
two observables which reveal complementary information on the mSUGRA parameter 
space: the "gold-plated" decays of charginos/neutralinos to trilepton final states and 
the rare decay B s — * ■ For the universal gaugino mass, M]/ 2 below 250 GeV 

and for the universal scalar fermion mass, Mq outside of 200-370 GeV we face a "no- 
lose" situation for the Tevatron: If tan/3 < 30, the Tevatron has a chance to see the 
trilepton events but not B s —* n + u~ during Run Ha, whereas for larger tan/3 the 
trilepton events become invisible (at least with an integrated luminosity of 2 fb _1 ) 
while B s — ► is enhanced to an observable level. For this study we perform an 

updated analysis of the trilepton signature which includes the full set of the decay 
matrix elements and spin correlations. This leads to a new, more promising search 
reach for the Tevatron. 



^ermilab is operated by URA under DOE contract No. DE-AC02-76CH03000. 
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1 Introduction and Motivation 

Even after imposing R-parity conservation by hand, the minimal supersymmetric standard 
model (MSSM) has an embarrassingly large number of parameters. The construction of 
a model with only a few parameters is desirable for aesthetic and practical reasons. The 
minimal supergravity (mSUGRA) scenario is currently the most elegant solution to 
this problem. In mSUGRA, local supersymmetry (supergravity) is spontaneously broken 
in the hidden-sector: the minimum of the scalar potential violates supersymmetry. The 
superpartner of the graviton, the spin 3/2 gravitino, acquires a mass, via this super-Higgs 
effect ||. Supersymmetry breaking is communicated to the observable sector via gravity 
interactions. The quantum gravitational corrections are in principle non-renormalizable, 
but taking the "flat" limit where the Planck scale goes to infinity (M P = (8irG N )- 1/2 -> oo) 
we are left with an effective theory which is renormalizable. The most general form of the 
low energy scalar potential after supersymmetry breaking contains the following free 
parameters:^ 

1. mass terms for the scalar particles, m^, 

2. trilinear couplings, hijk, 

3. the bilinear terms, fry. 

These terms break supersymmetry softly in the sense that no quadratic divergences appear 
in the theory The number of parameters is dramatically reduced if the superpotential 
splits into a hidden-sector and a visible sector part ||. Then these parameters are all 
universal 

rriij = M Sij, h^k = A Yij k , bij = B 8^. (1) 

Bq and the magnitude of the superpotential Higgs mixing parameter \x are fixed by elec- 
troweak symmetry breaking M. Outside of the scalar sector, there is a further universal 
soft supersymmetry breaking parameter: the common gaugino mass, Mi/ 2 - In the Higgs 
sector there is an additional parameter: tan /3, the ratio of the vacuum expectation values 
of the two neutral, CP-even Higgs fields. We thus have a total of 5 parameters in mSUGRA 

M , A , Mi/a, sgn/i, tan/3. (2) 

The parameters M , A , Mi/ 2 , are expected to be universal at the Planck scale, where 
supersymmetry breaking is mediated from the hidden-sector to the observable sector. In 
order to obtain weak-scale predictions we must use the renormalization group equations 
(RGEs) to evolve the parameters from Mp down to Myy. In addition, there could be 
unification at the grand unified (GUT) scale Mx = O(10 16 GeV) in an unknown simple 
group such as SU(5) or 5*0(10). This would strongly modify the RGEs between the 
Planck scale and the GUT scale. In order to avoid this model dependence, we consider 
M , A , Mi/2, to be universal at the GUT scale. We thus neglect any effects of the running 

2 If R-parity is conserved there is only one bilinear parameter = b. 
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of the RGEs between Mp and Mx- In our numerical analysis, tan (3 and sgn /i are fixed at 
the weak scale. 

The mSUGRA model with universal parameters is a simple, well motivated model for 
the low-energy supersymmetric spectrum. If an observable in this model deviates substan- 
tially from the standard model (SM), then we would expect to find a similar deviation in 
a more general non-universal case, modulo accidental cancellations. Thus we expect to be 
able to obtain more general results already in this simplified model. However, correlations 
between observables in mSUGRA such as B s — > versus (g — 2) M are typically lost 

in the more general case due to the additional parameter freedom. 

After more than 25 years there is still no experimental evidence for supersymmetry 
(SUSY). In principle we expect two kinds of signatures: 1) an indirect signature via a 
deviation from a SM prediction or 2) a direct detection of a supersymmetric particle. It is 
the purpose of this letter to demonstrate the complementarity of these two approaches in 
two specific examples. For the direct search we consider the trilepton signature, which was 
first proposed in the context of supersymmetry in [P-|Tl|, and has since become a standard 
search mode for supersymmetry production. We wish to investigate the correlation of this 
signature in mSUGRA with the rare B s meson decay B s — > fi + fi~ . We expect these two 
signatures to be observable at Run II of the Tevatron. 

In the following we present the actual status of leptonic B meson decays. We then 
proceed and perform a brand new analysis of the multilepton signatures in mSUGRA. 



We confirm the results of Matchev and Pierce [12|,[T3|]. We then go beyond this work to 
include the full matrix elements for the supersymmetric decays, which were not included 
in ISASUSY at the time. This leads to a more promising discovery reach for Run II of the 
Tevatron. 



2 The decay B s — > 

In the SM, the leptonic decays B s ^ — > £ + £~ are dominated by Z-boson penguin and box 
diagrams involving top quark exchange. The decay amplitude is helicity suppressed and 
thereby proportional to the lepton mass. The resulting small branching ratios are well 
below the upper bounds set by past and present experiments; the current situation is 
summarized in Table [I]. The uncertainty in the SM predictions |14| of Table [I] mainly 
stems from the decay constants f Bd = (200 ± 30) MeV and f Bs = (230 ± 30) MeV . 



While the element \V ts \ = 0.040±0.002 of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
entering the B s decays is well-known from CKM unitarity, the B s decays suffer from an 
additional uncertainty due to \Vu\- From a global fit to the unitarity triangle one finds 
0.15 < \V td /V ts \j^0.23 p[. We use the MS b-quark mass m h {m b ) = 4.25 GeV and 
assume a fixed^ MS value for the top quark mass, m t (m t ) = 167 GeV. This corresponds 
to a pole mass of m t = 175 GeV. There is an additional small uncertainty of ±0.3 x 10 -9 



3 In the mSUGRA analysis below, m t (m t ) is not fixed due to the variation of the strong coupling 
constant a s with the SUSY spectrum under the assumption of gauge coupling unification. However, the 
effect of this variation is small compared to the hadronic uncertainty. 
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Channel 


Expt. 


Bound (90% CL) 


SM prediction 14J 


B s -> e+e- 


L3 [17| 


< 5.4 x 10~ 5 


(8.9 ±2.3) x 10- 14 


B s -> 


CDF [TB| 


< 2.0 x 10~ 6 


(3.8 ± 1.0) x 10~ 9 


5 S -> r+r" 


LEP \V5\ 


< 0.05 


(8.2 ±2.1) x 10~ 7 


fid -> e+e" 


CLEO [0 


< 8.3 x 10~ 7 


(2.4 ±0.7 ±0.7) x 10" 15 




CLEO 


< 6.1 x 10~ 7 


(1.0 ±0.3 ±0.3) x 10" 1U 


^ -> r+r- 


LEP [19| 


< 0.015 


(2.1 ±0.6 ±0.6) x 10" 8 



Table 1: The experimental status and the SM predictions for the branching ratios B(B S ^ — > 
£ + £~). The error in the B s branching ratios mainly originates from the uncertainty in fs a , 
the two errors in the B d branching ratios correspond to the uncertainties in fs d and \V t d\- 

in B(B S — > fi + fi~ ) when scanning the top mass in its experimental range. Since we are 
interested in order-of-magnitude effects from supersymmetry, we can safely neglect the 
hadronic uncertainty and work with the central values of the decay constants and CKM 
elements. 

As B{B — > £ + £r) oc m|, the branching ratio is largest for £ = t. Yet r-lepton searches 
are very difficult at hadron colliders. Therefore the best mode to search for new physics 
is B s — > fi + fi~ . Still the decay to t + t~ has a very weak bound and we believe that an 
analysis of the existing LEP data could result in a stronger bound. Leptonic branching 
ratios of B d mesons are smaller by a factor of |V^/V^| 2 < 0.05 than the leptonic decays 
of B s mesons. From Table [l] we see that CLEO |20| has provided strong bounds for the 
electron and muon final states. Again, a dedicated search for the final state r's remains to 
be performed. The BaBar and Belle experiments could improve on Table [l] as well. 

From now on we restrict ourselves to the decay mode B s — > . In the SM, the 

decay mode B s — ► is experimentally challenging due to its small branching ratio. 

However, it has a quite distinctive and unique signature. During Run I of the Tevatron the 
CDF experiment has set an upper bound at 95% CL. [[L8[] 

B(B S -> n+ir) < 2.6 x 10~ 6 , (3) 

based on an analysis using very tight cuts to achieve an essentially background-free frame- 
work. When extrapolated to Run Ha, which corresponds to an integrated luminosity of 2 
fb _1 , this method yields a single event sensitivity of B(B S — > /U + /i - ) = 1.0 • 10~ 8 ||21|| . A 
recent analysis jZ^] found that CDF can discover B s — ► in Run lib with an integrated 

luminosity of 15 fb _1 , if B(B S — > /U + /i~) > 1.2 x 10 -8 . Currently further sophisticated stud- 
ies are underway, and better cuts may well improve the explorable range of B(B S —* fi + fi~). 
The SM range will certainly be accessible to the LHC experiments [^3 . 



The helicity suppression gives the leptonic B decays a very special role in the search for 
new physics, because they are sensitive to new chirality-flipping interactions like extended 
Higgs sectors with enhanced Yukawa couplings. This scenario occurs in the two-Higgs- 
doublet model (2HDM) and the MSSM, if tan (3 is large. In particular no other observable 
in B physics is more sensitive to effects from non-standard neutral Higgs bosons than 
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BiB — > £ + £~). In the type-II 2HDM the extra contribution to the decay amplitude grows 
like tan 2 /3 and is roughly of the same size as the SM amplitude [p4],|25||. In the MSSM, 
however, the amplitude involves three powers of tan /3, so that 

B(B S -+fi + iT) octan 6 /?. (4) 

This leads to an enhancement over the SM value by up to three orders of magnitude. This 
holds even in the flavour-universal MSSM, where flavour changes occur in weak vertices only 



and come with the usual CKM elements [25-27]. Moreover, in the limit of infinite SUSY 



mass parameters the MSSM contribution does not vanish, because the MSSM approaches 
a general 2HDM with flavour- changing neutral Higgs couplings in this limit. This feature 
and the resulting importance of B physics to probe the large tan/3 regime of the MSSM 



was first emphasized in [28]. The supersymmetric contributions enhancing B(B S — > ) 



will also affect the branching ratio and forward-backward asymmetry in the decay 

B -> (<ir,K)£ + r, (5) 

but these effects are small and strongly correlated with B(B S — > and will only be 

measurable, if B(B S — > is close to its experimental upper bound [p7|,|29|. 



The first analysis of B s — > in the minimal supergravity scenario (mSUGRA) 

revealed an interesting correlation of B(B S — > with the muon anomalous magnetic 

moment (g — 2)^ and the mass of the lightest Higgs boson H: If (g — 2)^ exceeds the SM 
prediction by 40 x 10~ 10 (which is the maximum allowed from BNL [[3(],|31| ) , B(B S — > 
is larger by a factor of 10-100 than in the SM and within reach of Run II of the Tevatron. If 
the decay B s — > fi + fi~ is observed at Run II of the Tevatron, then in mSUGRA the mass of 
the lightest neutral CP-even Higgs boson should be less than 120 GeV. The results of Ref. [[?]] 



were recently confirmed in Ref. [22], where also constraints from dark matter were included. 
A nice observation was made in Ref. |j32] : alternative SUSY breaking scenarios imprint 
signatures on (g — 2)^ and B(B S — > // + which are very different from those found for 
mSUGRA in Ref. 0. A substantial enhancement of B(B S — > in anomaly-mediated 

SUSY breaking scenarios is impossible and also hard to accommodate with gauge-mediated 
SUSY-breaking. Thus an observation of B s — > fi + fi~ at Run II of the Tevatron can be 
viewed as a signature of supersymmetry with gravity-mediated supersymmetry breaking. 



In our numerical analysis of B(B S — > fi + fi ) we use the formulae of [^7| supplemented with 



an all-order resummation of tan /5-enhanced supersymmetric QCD corrections p3|. Other 



QCD corrections (neglecting gluino contributions) were calculated in |3^] and found to 
be small, if the renormalization scale is chosen close to the masses of the SUSY particles 
propagating in the loop. 



3 Multilepton SUSY events 



Trileptons are an old phenomenological signature of cascade decays in particle physics 
|35|,|36||. They were first applied in supersymmetry in the cascade decays of electroweak 
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gauge bosons via light gauginos |||3 A first systematic study was performed in P-[TT 
It was first shown in |38] that off-shell W, Z cascade decays via gauginos can also lead 
to promising trilepton signatures. The study was however restricted to tan/? = 1. Also 
the sparticle spectrum was not computed in terms of a small set of MSSM parameters. 
These points were generalized in |39|]. An extension of the Tevatron studies to the LHC 
was performed in pE0|]f| and to a yfs 



4TeV Tevatron in W2 



luminosity Run II at the Tevatron was performed in |1-5|.| 1 1 
call mSUGRA was performed in f 



A first study for the high- 
A first study of what we now 
including a multichannel analysis. Various unification 

In 13 ~ 



scenarios were probed via the trilepton signature in |46|,|47|]. In ||| the case of large 
tan /3 was first considered. This can lead to additional contribution of r leptons in the 
final state which was discussed in detail in [p0| -|5^]. An update improving many features 
of the analysis has recently been performed in Hl2, 55fl including a reevaluation of the W^* 
background fl^fl. We now discuss in detail how we have improved on this work. 



3.1 Numerical Analysis 



In general our strategy for the analysis of the multilepton signatures closely follows that of 
Matchev and Pierce |T2],|I1|. The major differences between these studies are: 



While Matchev and Pierce [0,0] used leading-order (in QCD) cross sections for 
both the signal and background processes we use next-to-leading order calculations 
for the background and most important signal processes. 



In Refs. [|T2| , p~3j| the authors used ISAJET 7.42 |56| for the simulation of the signal. For 
the background they used PYTHIA 6.115 |57|| . In Ref. JT3| they used a combination of 
PYTHIA and COMHEP @ for the WZ background.^ Instead, we use HERWIG6.4 
[j9H to simulate the signal and a combination of HERWIG and COMPHEP to simulate 
the background. This has a number of advantages due to recent improvements in 
HERWIG. While ISAJET7.42 uses matrix elements for the three body decays of the 
gauginos it performs the production and decay of the SUSY particles independently. 



HERWIG6.4 makes use of the method described in |6(J to perform the SUSY decays 
including all the momentum and spin correlations. One of the major advantages 
of our studies is that the spin correlation algorithm used in HERWIG can use the 
TAUOLA decay package pT| to give the correct helicity of the tau on an event-by- 
event basis, rather than averaging as was done in |T2,[r^]. ISAJET7.58 is used to 
calculate the SUSY spectrum and decay rates. 



In our simulation we make use of the PGS detector simulation. This is a more recent 
version of the SHW package. The default PGS parameters are used together with the 
following additions. 



4 Implicitly these signatures appear in Ref. |37[| . 

5 For analogous signatures for gluino production see for example |4l|| 

6n 



3 This brief summary is based on a private communication with Konstantin Matchev and differs slightly 
from the presentation in the papers. 
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1. A requirement that the total transverse energy in a cone of size AR = ^A<p 2 + Ar] 2 = 
0.4 about the direction of a muon is less than 2 GeV. Here is the azimuthal angle 
of the particle, r\ = — lntan(#/2) is the pseudo-rapidity, and 6 is the angle of the 
particle with respect to the incoming beam. 

2. The missing transverse energy, Et, which is calculated in PGS using the calorimeter 
cells, is corrected to include any muons. The transverse energy E? = Esiri8, where 
E is the energy of a particle. 

3. A cut on the ratio of the electromagnetic, E em , to hadronic energy, -Ehad, in a jet 
E em / Eh a d < 10 to avoid electrons being considered as jets. 



These modifications are virtually identical^] to those made in Ref. [13]. 



3.2 Background and Signal 

The major backgrounds to multilepton production are: 

WW : Here the background is from the production of leptons in the cascade decay of the 
W bosons. These events are simulated using HERWIG6.4 and the results normalized 



to the next-to- leading order result from MCFM [63] 



WZ : In this case the leptons come from the decay of the gauge bosons. This is the main 
physics background in most of the channels. As this channel is so important, rather 
than use HERWIG, which does not include all the diagrams, we use COMPHEP to 
generate all the Feynman diagrams for the four-fermion final state. HERWIG is then 
used to generate the initial-state radiation and hadronization of the partonic events 
generated with COMPHEP. The results are normalized to the next-to-leading order 
cross section calculated using MCFM |)3[ . 



ZZ : The background comes from the production of leptons in the decay of the Z bosons. 
These events are simulated using HERWIG6.4 with a next-to-leading order normal- 



ization from MCFM 63 



W+jets: In these events the background comes from one real lepton produced in the decay 
of the W and additional fake leptons from the mis-identification of isolated hadrons 
as leptons. The simulation of isolated hadrons producing fake electrons and muons 
cannot be performed reliably using Monte Carlo event generators and therefore we 



adopt the same procedure as in Ref. |L3[ in order to estimate this rate from the 



Run I CDF data. The SHW simulation, which was the predecessor to PGS, gives 
a reasonable simulation of fake tau production []13"|, p4}] , and we therefore rely on the 
Monte Carlo simulation for fake tau production. HERWIG6.4, which includes the 



matrix-element correction to the parton shower as described in Ref. [65|, is used to 



7 In Ref. |f3| the tracking coverage is extended from \r]\ < 1.5 to \rj\ < 2.0 however this is now the default 
in PGS. 



s 



simulate this background. The results of the Monte Carlo simulation are re-scaled to 



the next-to- leading-order cross section calculated using MCFM ]66 ] . 



Z+jets: In this case the background comes from either the leptons produced in the cas- 
cade decay of the Z, with additional fake leptons as described above, charge mis- 
identification, or leptons from the Z itself depending on the channel. Again HER- 
WIG6.4 which includes the matrix-element correction to the initial-state radiation is 
used to simulate this process. As with the W+jets background, the results of the sim- 
ulation are normalized using the next-to-leading-order cross section calculated with 



MCFM ||66|| . The fake lepton rate is calculated as for the W+jets events. 



tt : The background comes from the semi-leptonic decay of the top quarks. We simulate 
this background using HERWIG6.4 and use the next-to-leading-order with next-to- 



leading-log resummation calculation of Ref. [67] to normalize the results. 



There is also a potential background from the production of bb followed by the semi-leptonic 
decay of the produced b-flavoured hadrons. However, this should be negligible after the 
isolation and cuts on the transverse momentum, p?, of the leptons. 

In all cases we use HERWIG6.4 to simulate the signal. The cross sections for xtxt an d 
XiXi) which give the dominant contribution to the SUSY production, are normalized to 
the next-to-leading order cross section using the results of [|6^] . The leading-order cross 
section from HERWIG is used for the remaining processes. 



3.3 Cuts 

The basic idea of the analysis is to take several different cuts, which should differentiate 
between the signal and the background, consider all possible combinations of the cuts at 
each point in parameter space and optimize the choice of cuts in order to maximize S/ yB, 
where S and B are the number of signal and background events respectively. In addition 
to the requirement of S/ \f~B~ > 5 all the contours shown also require that there at least five 
signal events. 

The cuts can be split into two types: channel-independent cuts and those which are 
specific to a given channel. 

The channel independent cuts we use are: 

1. Five possible cuts on the missing transverse energy, > {15, 20, 25, 30} GeV or no 
cut; 

2. An optional veto on the presence of QCD jets with Et > {10, 15, 20, 25, 30} GeV in 
the event. 

3. A cut on the minimum invariant mass of any opposite sign same flavour electron or 
muon (OSSF) pairs in the event m u > {5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60} GeV, 
or no cut. 
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4. A cut on the maximum invariant mass of any OSSF pairs in the event 
rriu < {50, 60, 70, 80} GeV, or a cut on the minimum difference between the mass 
of any OSSF pair and the Z mass |M Z — m u \ > {10, 15} GeV, or no cut. 

5. A cut removing any events which have an electron or muon with trans- 
verse mass 65 GeV < hit < 85 GeV. The transverse mass is given by 
m% = 2\pTe\ |pTmiss| (1 — cos A0), where pti is the transverse momentum of the lepton, 
PTmiss is the missing transverse momentum and A0 is the azimuthal angle between 
the lepton and the missing transverse momentum. 



The channel dependent cuts are: 



1. For the trilepton and like-sign lepton signature we require a central electron or muon 
with 1 77 1 < 1 and pt > 11, 15, 20 GeV, or no requirement for a central lepton. 

2. For the lepton+ tau jet channel we require that the tau decay only has one prongP] and 
that the fraction of the tau jet's momentum carried by the charged track is greater 
than 80%, as suggested in |69| . 



3. A cut on the transverse momenta of the leptons is also applied depending on the 
channel, the cuts for the various different channels are given in Table 0. 

These cuts are the same as used in Ref. |12|,[13| apart from one extra missing transverse 
energy cut and the cut on the momentum of the track in the tau jets. 



4 Discussion 



4.1 Trilepton (3L) Channel 

We first discuss the trilepton (3L) channel. The production cross section for the chargino- 
neutralino pair, xfx2°, scales like M x j^ 2 and is thus only relevant in regions with small 

M1/2, i- e - M1/2 ~ 300 GeV. We present our results in Fig. 1 on the M — Mi/ 2 plane, for 
two different values of tan/3: 5 and 50. Starting with the (3L) events in the small tan/3 
plot, we observe that the Tevatron could reach a 5<r discovery depending on its luminosity 
in the regions 



Region I (3L, tan (3 — 5) 



Region II (3L, tan f3 = 5) 



3 i.e. the tau jet only had one charged track. 



M 


< 


175 GeV 


, M1/2 


< 


230 GeV 




2/6" 


1 


M 


< 


190 GeV , 


M1/2 


< 


270 GeV , 


£ = 


10/6- 


-1 


M 


< 


200 GeV , 


M1/2 


< 


290 GeV , 


C = 


30/6- 


-1 


M 


> 


480 GeV , 


M1/2 


< 


190 GeV , 


C = 


10/6- 


-1 


M 


> 


360 GeV , 


M1/2 


< 


220 GeV , 


£ = 


30/6- 


-1 
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Trilepton 


channel 




Cut Set 




P T {k) 


PT(h) 


1 


11 


5 


5 


2 


11 


7 


5 


3 


11 


7 


7 


4 


11 


11 


11 


5 


20 


15 


10 


Like-sign dilepton channel 


Cut Set 


Pt(£i) 


p T {h) 




1 


11 


9 




2 


11 


11 




3 


13 


13 




4 


15 


15 




5 


20 


20 




Dilepton+hadronic tau c 


lannel 


Cut Set 


Mk) 


p T (h) 


Pt(t) 


1 


8 


5 


10 


2 


8 


5 


15 


3 


11 


5 


10 


4 


11 


5 


15 


5 


15 


10 


15 



Table 2: Lepton pt cuts, all the cuts are given in GeV. 



Region I is larger than in Ref. |~2],[T3| due to the extra missing transverse energy cut we 



have made. Region II is larger than in W\ |T^] due to the lighter spectrum produced by 
ISAJET7.58 which gives a larger signal cross section than for the same SUGRA point with 
ISAJET7.42. Overall, we thus obtain a more promising trilepton search reach at Run II of 
the Tevatron than in previous studies. This is one of the main conclusions of this paper. 
Although we present our analysis in terms of the fundamental parameters of mSUGRA 



Mq, Mi/2, tan/?, there exist mass sum rules |7D[ relating these parameters to the physical 



masses. We can thus approximately translate the mSUGRA parameters into the physical 
masses (up to a 10% accuracy), 

2m^o ~ m^o ~ m~± ~ -m„ ~ 0.8Mi/ 2 , 



m 5 ~ ^M 2 + 4M 1 2 /2 , (6) 

where m^o are the neutralino masses with the lightest being the LSP. m-± denote the 
chargino masses. m g and vtiq denote the gluino and the heavy first and second generation 
squark masses, respectively. Thus the Tevatron trilepton reach is relevant for a lightest 
chargino and a next-to-lightest neutralino mass below 200 GeV and (indirectly) for a gluino 
mass below 600 GeV. 
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tan/S=5, A = 0, / a>0, m t =l75 GcV tan/3 = 50, Aq = 0, /zX), m t =l75 GeV 




200 400 600 800 1000 200 400 600 800 1000 



M [GeV] Mq [GeV] 

Figure 1: Tevatron 5a reach in the trilepton (3L) channel in the M — Mi/ 2 plane in the 
mSUGRA scenario for small tan/? = 5 (left) and large tan/3 = 50 (right) and A = 0, 
fi > and m t = 175 GeV. The shaded areas indicate the 5a reach for an integrated 
luminosity of 30 fb^ 1 (magenta), Wfb^ 1 (blue) and 2fb~ l (green) and from top to bottom, 
respectively. Dashed contours represent the SUSY contribution to the muon anomalous 
magnetic moment (in units of 10~ 10 ) and the dotted contours are iso-mass contours of the 
lightest neutral Higgs boson. The solid contour ( only for tan (3 = 50) indicates the prediction 
for the branching ratio B{B S — > fi + /j~). In the left (tan/3 = 5) plot the solid line indicates 
where itl^± = and the dot-dashed line itl^± = . The large red regions are excluded 
by theory and experiment /Q/. 



In Fig. 1 we include Higgs boson mass isocurvatures. The calculation of the light Higgs 
boson mass strictly follows the procedure described in Ref. fH|. In Region I the light Higgs 



boson mass varies from (106 4- 112) ± 3 GeV where the error denotes a theoretical error in 



the calculation |[72|| . The current status for the Higgs boson mass bound from LEP is [73 
Mh > 113.5 GeV for sin 2 (/3 — a) ~ 1, where a is the mixing angle in the CP-even Higgs 
sector. Thus tan (3 = 5 has to be considered as the minimum tan (3 value where trilepton 
events can appear at the Tevatron since lower values will be in conflict with the LEP Higgs 
bound. The Higgs mass increases with tan f3. 

For small values of tan/? the super symmetric contribution to the muon anomalous 
magnetic moment (<5a^ USY ) is in general small [ 74| . However, observing a trilepton signal 



in Region I will imply an at least moderate enhancement of a M by 

Saf SY > 7 x 10- 10 , (7) 

due to the light charginos and sneutrinos. This is not the case in Region II, where £a^ USY 
approaches zero. 

There is a "null" M — M\i2 region at the Tevatron where no (3L) signal can be observed. 
This "null" region spans the parameter space where 200 < M < 370 GeV and is due to 
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the fact that off-shell, slepton mediated decays of the gauginos destructively interfere with 
the gauge boson mediated decays. Since the decay B s — > ji + cannot be seen for such 
small values of tan (3 the "null" area has to be covered by other Tevatron searches (possibly 
Higgs searches J75) ). 

In Region II, a ha signal reach can be achieved with C = 10,30fb _1 . The light Higgs 
boson mass varies between (106 -j- 111) ± 3 GeV and the supersymmetric contributions to 
the muon anomalous magnetic moment are small, 

6af SY < 5 x 10- 10 . (8) 

The CLEO/Belle/ALEPH allowed band for the branching ratio B(b — * 57) does not 
set any constraint in the low tan (3 = 5 region. In what follows, we adopt the conservative 
value®] 2 x 10" 4 < Bib -> 57) < 5 x 10" 4 . 

We now turn to the discussion of the large tan/3 = 50 plot in Fig. 1. The new Region I 
(low Mo) is not there any more (see below). The new Region II is given by 

n • tt /or x „ / M > 480 GeV , M 1/2 < 185 GeV , £ = 10 fb' 1 , 

Region II (3L,tan/3 = 50) : ( Mq > 36Q GeV Mi/ 2 < 230 GeV , £ = 30/^. 

The trilepton search only covers a small range of the Mo — M1/2 plane. There is a further 
strong constraint. The light Higgs boson mass varies in the range (110 4- 119) ± 3 GeV 
from bottom to top where the contour = 115 GeV is indicated in Fig. 1 as a dotted 
line. The LEP preliminary bound, > 113 GeV, is evaded (almost) everywhere in this 
plot. However, this is not true for the b — > 57 constraint. Using the results in Ref. [75] we 
find that the excluded region due to the b — > S7 constraint (accidentally) coincides almost 
exactly with the area below the = 115 GeV Higgs mass contour. We thus do not redraw 
it. This leaves only a small area for the 3L events to appear at the Tevatron 

Region II' (3L, tan/3 = 50) : M > 800 GeV ; M 1/2 ~ 190 ^ 230 GeV, C = 30/r 1 . 

Notice also that this area is bounded for large Mo due to the radiative electroweak symmetry 
breaking requirement. 

By contrast, the decay mode B s — > fi + ^~ is very promising for tan/3 = 50. The 
branching ratio varies from a handful of events for large M up to 100 events for low M 
already with L = 2fb~ 1 . This very nicely complements the trilepton search and is the 
second main result of our paper: the tri-lepton search covers most of the low tan (3 region, 
the branching ratio B{B S — > /i + /U~) covers most of the high tan/3 region. 

Fig.l also verifies what had been seen in Ref. M: the branching B(B S — > /i + /U~) and 
the anomalous magnetic moment of the muon, <5a^ USY , go hand-in-hand, i.e. there is 
a strong correlation. New experimental results for the anomalous magnetic moment of 
the muon are expected to appear as this manuscript is written. The up-to-date BNL 
excess is -4.2 x 10~ 10 < 5(1^ < 41.3 x 10~ 10 , at 90% C.L |J{J[TT| with the central value 
So,™ SY = 18 x 10- 10 0. 



We follow here the discussion in the footnote no. 18 of Ref. 176 
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Figure 2: Tevatron 5a reach in the trilepton (3L) channel in the M\m — tan/5 plane for the 
mSUGRA scenario with a) M = 150 GeV, b) M = 300 GeV and c) M = 600 GeV and 
A = 0, fi > and m t = 175 GeV. The shaded areas indicate integrated luminosities of 
30fb~ x (magenta), Wfb^ 1 (blue) and 2fb~ l (green) and from top to bottom, respectively. 
Dashed contours represent the SUSY contribution to the muon anomalous magnetic moment 
(in units of 10 and the dotted contours the Higgs boson mass. The solid contours 
indicate the prediction for B(B S — ► 



Our next step is to present the simulation of the 3L events in the Mx/2 — tan/3 plane 
for different values of M . Our results are depicted in Fig. 2. For M = 150 GeV, the 5cr 
reach of the Tevatron with C = 30/6 -1 is bounded by tan/? < 30 and Mi/ 2 J$ 270 GeV. 
In this region the Higgs mass is always below 114 GeV. With an integrated luminosity of 
£ = 2fb~ 1 , the Tevatron can observe the 3L events in a region with small tan/3 < 5 and 
Mi/2 < 210 GeV. In this region, the maximum Higgs mass is 109 GeV. 
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tan/S=5, A = 0, fx>0, m t =l75 GcV tan/5=5, A Q = 0, /x>0, m t =i75 GcV 




M [GeV] Mq [GeV] 



Figure 3: Tevatron reach in the like sign dilepton channel (2L) (left) and and dilepton plus 
tau jet channel (2L1T) (right) for low tan/3 = 5 and all the other parameters as in Fig.l. 
We do not observe these two modes at large tan/3 = 50. Contours and shaded areas are as 
in Fig.l. 



As tan (3 increases, the stau mass decreases and the gaugino decays are mediated mostly 
by staus rather than sleptons. Gradually the 3L event rate is decreasing as we increase 
tan/3. When the 3L events are no longer observable, the branching ratio B(B S — > 
becomes significant and the decay B s — > has an observable rate of up to 30 events with 

£ = 2fb~ 1 . Moving to higher values of M = 300 GeV (see Fig.2), the 3L events dissappear 
for the reason we have already explained when discussing Fig.l. However, a large portion 
of the parameter space with tan /3 > 30 is still accessible to B s — > fi + fi~ searches. The 
parameter space with tan/? values below 30 has to be covered by other searches, most 
probably Higgs searches (squarks and gluinos are too heavy to be directly discovered at 
Tevatron in the mSUGRA scenario |77|.) The two cases M = 150 GeV and M = 300 GeV 
have to be considered respectively as the most "promising" and the "nightmare" benchmark 
scenarios for the Tevatron Run II. In order to complete the discussion on M1/2 — tan (3 plane 
we show in Fig. 2c the corresponding plot for Mq = 600 GeV. This region is intermediate 
in search sensitivity and is bounded by Mi/2 < 210 GeV. It is largely independent of 
tan /3. A considerable amount of the remaining parameter space is covered by the search 
from B s — > fi + fi~ . We observe a considerable overlap between the two search regions for 
tan /3 > 40. For a complete analysis of the B s — > fi + {i~ mode and other observables on 
M1/2 — tan (3 plane the reader is referred to Refs. 171 p2 . 



4.2 Other than 3L Channels 



Apart from the trilepton (3L) analysis, we also study the modes with like sign dileptons 
(2L), dilepton plus tau-jet (2L1T), lepton plus two tau-jets (1L2T), three tau-jets (3T), like 
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sign lepton plus tau-jet (1L1T), and two like sign tau-jets (2T) in the final state. The idea 
is to identify alternative signatures in regions of the parameter space where the 3L signal 
is suppressed. For all these channels only the 2L and 2L1T survive our cuts and only in 
the small tan/5 ~ 5 region. The results are shown in Fig. 3. Depending on the luminosity, 
we observe a 5a reach at the Tevatron, in the mSUGRA regions: 



r M < 120 GeV , M va £ 175 GeV , £ = 2/V 1 , 
(2L,tan/3 = 5) : I M < 150 GeV , M 1/2 < 200 GeV , £ = 10/Tr 1 , 
^ M < 160 GeV , M 1/2 < 230 GeV , £ = 30/r 1 , 

f 40 < M < 80 GeV , M 1/2 < 180 GeV , £ = 2fb^ , 
(2LlT,tan/3 = 5) : \ 50 < M < 100 GeV , M V2 < 260 GeV , C = 10 fb' 1 , 

^ M < 120GeV , M 1/2 < 230 GeV , £ = 30/6" 1 , 

Comparing with the Region I of the 3L signal, we see that an observation of 2L or/and 
2L1T will pin down the parameters Mo — M\i 2 . In the Mi/ 2 — tan/3 plane and for the 
optimistic benchmark M = 150 GeV, 2L events are only observable with £ = 10(30)/6 _1 
for Mi/ 2 < 190, tan (3 < 3 (M 1/2 < 220, tan f3 < 6). 

The search reach we obtain for signatures involving more than one hadronic tau are 
worse than those obtained previously |H| . The main reason for this is that we started from 
the leading-order Drell-Yan production processes with additional jets generated by the 



parton-shower algorithm for W and Z production, whereas [69] started from the W/Z+jets 



process. The study in |69[] was also at the parton-level rather than using fast detector 
simulation. This means we have a much larger background from W/Z with fake taus. This 
agrees with the results of []53] where the signals with at least two electrons/muons gave the 



largest reach. It should also be noted the we took no account of the problems of triggering 
on the multi-hadronic tau final states which will make the situation even worse. 



5 Conclusions 



We have presented two complementary observables to study the parameter space of minimal 
supergravity at Run-II of the Fermilab Tevatron: the production of charginos or neutralinos 
decaying into trilepton (or other than trilepton) final states and the branching ratio of the 
rare decay B s — > fi + fi~ . The simulation of the trilepton events has been improved over 
previous results in various respects. In particular next-to-leading QCD corrections have 
been incorporated where available and the full momentum and spin correlations between 
production and decay of the supersymmetric particles have been taken into account. We 
have also presented modified cuts. With these enhancements we find the Tevatron reach for 
trilepton events increased. The two studied observables reveal very useful complementary 
information on the mSUGRA parameter space, covering the regions of small and large 
tan/3, respectively. Our results can be read directly from the Figs. 1,2. For example, for 
Mo = 150 GeV, Aq = and fi > we find that a Run-IIb with 30 fh~ x will probe values of 
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Mi/2 in the range 

tan 3 

M 1/2 < 250 GeV - 70 GeV x — 

which decreases with tan/3. On the contrary, if the Tevatron can measure B(B S — ► fi + n~) 
down to 10~ 8 (which requires less than 30 of integrated luminosity), this will probe 
the region with tan 8 > 32 and the range 

M 1/2 < 19 GeV x tan B - 260 GeV 

for the chosen values of M , A and /i. Similar results can be drawn for other choices of 
the parameter Mq. 

This work should be considered as a first attempt of relating "direct" SUSY searches 
such as the 3L events with an "indirect" rare B-decay, such as B s — > at the Tevatron. 
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